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Abstract 

We consider an extension of the Standard Model in which the symmetry is en- 
larged by a global flavour factor A4 and the scalar sector accounts for three copies of 
the Standard Model Higgs, transforming as a triplet of A4. In this context, we study 
the most general scalar potential and its minima, performing for each of them a model 
independent analysis on the related phenomenology. We study the scalar spectrum, 
the new contributions to the oblique corrections, the decays of the Z and W^, the 
new sources of CP and flavour violation, which all are affected by the introduction 
of multiple Higgses transforming under ^4. We find that this model independent 
approach discriminates the different minima allowed by the scalar potential. 
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1 Introduction 



The current data on neutrino oscillations seem to point at one small and two large angles 
in the neutrino mixing matrix (l]-[6] . The data are consistent with various mixing patterns, 
where in particular the agreement with the tri-bimaximal (7||8] mixing pattern is striking 

The use of non-Abelian discrete flavour symmetries has been proposed in different 
models (for a review see [9j) to generate both the mentioned lepton mixing patterns and 
the quark ones. In general, in those models, one introduces so called fiavons, scalar fields 
charged in the flavour space, usually very heavy. Once the fiavons develop specific vacuum 
expectation values (vevs), this translates to structures in the masses and mixings of the 
fermions. However, imposing the correct symmetry breaking patterns on the fiavons is 
highly non-trivial. This holds in particular if two or more fiavons are used, breaking in 
different directions in flavour space. So far, only a few techniques have been developed, 
all of which need a supersymmetric context or the existence of extra dimensions 

Alternatively, one can look at models that require only one flavour symmetry break- 
ing direction. In this case the scalar potential that implements the breaking can be non 
supersymmetric and does not require extra dimensions. Of particular interest is the pos- 
sibility that one set of fields simultaneously takes the role of the fiavons and the Standard 
Model (SM) Higgs fields, identifying the breaking scales of the electroweak and the flavour 
symmetries. 

In this paper, we will consider the discrete flavour symmetry ^4 and we will assume that 
there are three copies of the Standard Model Higgs field, that transform among each other 



as a triplet of A4 10 - 16 . The presence of this extended Higgs sector has an deep impact 
on the high energy phenomenology: indeed new contributions to the oblique corrections 
as well as new sources of CP and flavour violation usually appear in this context. We will 
analyse the constraints coming from these observables for all the vacuum configurations 
allowed by the scalar potential and will discuss on the viability of each of them. 

The structure of the paper is as follows. In section |2j we will introduce the scalar 
potential invariant under and under the gauge group of the Standard Model. In section 
|3] we will introduce the various physical Higgs fields that are present in the model. 

In the subsequent two sections we will present the different minima allowed by the 
potential and discuss the corresponding Higgs spectrum. These minima may or may not 
violate CP symmetry, depending on the question whether all vacuum expectation values 
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of the Higgses are real or not. In section |4] we will discuss the cases that do not violate 
the CP symmetry, while in section |5] we discuss CP breaking minima. 

Section [6] we will discuss bounds on the allowed parameters using respectively unitarity 
constraints, decays of the Z and bosons and constraints by oblique corrections. We 
note that all these bounds are rather model independent, meaning that they depend on 
the flavour symmetry assignment of the relevant Higgs fields, but not on those of the 
fermions in the theory. Further bounds can be derived from fermion decays and meson 
oscillations, but these bounds are always model dependent. We will present some of these 
in an accompanying paper 
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Finally, in section [7] we present the results of our analysis and in section |8] we conclude. 
In the appendix A we report useful formulae for the analysis of the T, S and U parameters. 



2 The A4 Scalar Potential 

We consider the Standard Model gauge group S'[/(3)c x SU{2)l x f/(l)y with the addition 



of a global flavour symmetry A4 [18, 19 . We consider three copies $a, a = 1,2,3, of 
the conventional SM Higgs field (i.e. a singlet of SU{3)c, doublet of SU{2)l and with 
hypercharge Y = 1/2) such that the three Higgses are in a triplet of the flavour group A^. 
Once the flavour structure of the quarks and leptons is specified, each will couple to 
the three fermion families according to the group theory rules in a model dependent way. 



We will study these couplings in more detail in 17 . 

Below, we will write down the most general scalar potential for the three Higgses that is 
invariant under the flavour and gauge symmetries of the model. After the fields occupy one 
of the minima of the potential, electroweak symmetry gets broken (while electromagnetism 
is conserved) and we can develop the fields around their vacuum expectation values as 

^ ^ ^ /Re$i+2lm<l>A ^ ^ / Re0i + ilm</)i \ 
I Re $0 + 2 Im $° i V2 L^e*'^'" + Re 0° + i Im0° i ' 

Here Vae^^"" is the vacuum expectation value of the a*^ Higgs field. One or two of the f „ can 
be zero, implying that the corresponding Higgs field does not develop a vev. Furthermore, 
if all vevs are real (so if all Ua are zero) CP is conserved, while if one or more are 
nonzero, CP is broken. Note that in general, there is the freedom to put one of the phases 
to zero by a global rotation. 

We will use the ^4 basis as developed by Ma and Rajasekaran (MR) flO). The analysis 



could also be done in a different A4 basis, for instance the one of Altarelli and Feruglio [20] . 
The results would then look different, but would obviously be equivalent. In the MR basis, 



the most general potential ^[$a] can be written as 10,21 



As 



+ 



+ e 



(2) 



in agreement with the usual notation adopted in the two Higgs Doublet Models (2HDM). 
The parameter /x^ is typically negative in order to have a stable minimum away from 
the origin. All the other parameters, Aj, are real parameters which must undergo to 
the condition for a potential bounded from below: this forces Ai and the combination 
Ai + A3 + A4 + A5 cose to be positive. 

It is interesting to notice that, contrary to other multi Higgs (MH) scenarios, here we 
can not recover the SM limit, with one light scalar and all the others decoupled and very 
heavy. The flavour symmetry constrains the potential parameters in such a way that the 
scalar masses are never independent from each other. This can be easily understood by 
a parameter counting: the scalar potential in eq. ^ presents 6 independent parameters 
and the number of the physical quantities is 8, i.e. the electroweak (EW) vev and the 
seven masses for the massive scalar fields. 

We will study the minima of the potential in eq. ([2]) under electromagnetism conserving 
vevs as specified in eq. ([T]) by studying the first derivative system 



a$7 



0. 



(3) 



where $x is of the fields Re Re Im $^ or Im $° and by requiring that the Hessian 

^^^W (4) 

has non negative eigenvalues, or in other words that all the physical masses are positive 
except those ones corresponding to the Goldstone bosons (GBs) that vanish. 

In sections |4] and |5] we will verify that this potential presents a number of solutions. 
Some of them are natural in the sense that they do not require ad hoc values of the potential 
parameters; these are only constrained by requiring the boundness at infinity and the 
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positivity of all the physical scalar masses. The only potential parameter constrained is the 
bare mass term /i^ which is related to the physical Electroweak (EW) vev, = 
Others require specific relations between the adimensional scalar potential parameters and 
may have extra Goldstone bosons. 



3 The Physical Higgs Fields 



The symmetry breaking of the Higgs fields of equation eq. ([T| leads to a large number of 
charged and neutral Higgs bosons as well as the known Goldstone bosons of the Standard 
Model. 

In the most general case, where CP is not conserved, the neutral real and imaginary 
components of eq. ([T]) mix to five CP non-definite states and a GB: 



IT" 



t/6(a+3)Im0° 



(5) 



Here a 
reads 



1, 2, 3 and a = 1 — 5, while a = 6 represents the GB 7r°. In matrixform this 

/Re0o\ 



u 



Re 00 
Im 



yim0°y 



(6) 



Clearly eq. ([s]) holds also in the CP conserved case: in that case the 6 by 6 scalar mass 
matrix reduces to a block diagonal matrix with two 3 by 3 mass matrices leading to three 
CP even states and 2 CP odd states and the GB vr*^. 

The three charged scalars mix into two new charged massive states and a charged GB. 



(Ht\ 



(A 



S 



(7) 



where 7r+ is the Goldstone boson eaten by the gauge bosons . In general, the S" is a 
complex unitary matrix. In the special case where CP is conserved, its entries are real 
(and it is thus an orthogonal matrix). 
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4 CP Conserved Solutions 
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In this section, we will study minima of the potential in eq. ^ in which only Re0° 
develops a vev, i.e. the CP symmetry is conserved. In this case we expect having 3 
neutral scalar CP-even states, 2 CP-odd states and 2 charged scalars as well as a real and 
a complex GBs originating from respectively the CP-odd states and the charged states. 

In the CP-conserved case all the cOa vanish and the first derivative system in eq. (|3| 
reduces to 

vi[2{vl + vi + vl)\i + {v^ + v'l){\s + \i + Xr, cose) + 211"^] = 0. 
W2[2(t;2 + t;| + t;|)Ai + K + f|)(A3 + A4 + A5Cose) + 2/i2] = 

V3[2{vl + vi + vl)\i + {vf + vj){\3 + \4 + \5c0se) + 2fx'^] = 0. 

— tiDAs sin e = 0, 

'y2('i^i — "^1)^5 sin e = 0, 

^3(^^| - ^^?)A5sine = 0, 

where the first three derivatives refer to the real components $° and the second ones to 
the imaginary parts. In the most general case, when neither e nor A5 is zero, the last three 
equations allow two different solutions 

1) vi = V2 = V3 = V = Vyj/ a/S; 

2) 7^ and f 2 = f 3 = (and permutations of the indices); in this case vi = Vw 
Both these solutions are solutions of the first three equations as well, provided that 

fi^ = — (3Ai + A3 + A4 + A5 cose)f^/3 for the first case 

(9) 

/i^ = — Aif^ for the second case. 

In this cases A5 can be chosen positive, as a sign can be absorbed in a redefinition of e. 

Next, we consider the case where sine is 0. This implies e = or vr. We may however 
absorb the minus sign corresponding to the second case in a redefinition of A5 that is now 
allowed to span over both positive and negative values. 

Assuming vi 7^ 0, we may solve the first equation in eq. ([s]) with respect to /x^. Then 
by substituting /x^ in the other two equations we get 

(10) 

V3{vj-vl){Xs + XA + X5) = 0. 



Next to the two solutions present in the general case, this system has two further possible 
solutions 

3) = 0,V2 = vi = Vui/V2 and permutations. This requires 

^2 = _ (4A^ + A3 + A4 + As) vl/A . (11) 

4) (A3 + A4 + A5) = 0. This condition implies that in the real neutral direction there is 
a enlarged-0(3) accidental symmetry that is spontaneously broken by the vacuum 
configuration, thus we xpect extra GBs. Indeed in this case vi, V2 and are only 
restricted to satisfy fi + f | + f 3 = and the parameter /i^ is given by /i^ = — Ai^;^. 

Finally, the case A5 = allows special cases of the solutions 1) to 4), but does not give 
rise to new solutions. For this reason, we will discuss only the general cases and the case 
e = in the remainder of this section and comment what happens for A5 = 0. 

4.1 e 7^ 0: The Alignment {v^v^v) 

In the basis chosen, the vacuum alignment (f , f , v) preserves the Z3 subgroup of It is 
convenient to perform a basis transformation into the Z3 eigenstate basis, 1, 1' ~ cu, 1" ~ cu^ 
according to 

^ = ($1 + $2 + <I>3)/V3 ~ 1 

ip' = ($1 + u;<I>2 + w^$3)/V3 ~ w 

ip" = ($i + a;2$2 + a;$3)/V3~a;^ (12) 

When ^4 is broken to Z3 in the Z-^ eigenstate basis, ~ 1 behaves like the standard Higgs 

doublets: its neutral real component develops a vacuum expectation values (v?^'^) = Vw 

and all its other components correspond to the GBs eaten by the corresponding gauge 

bosons. The physical real scalar gets a mass given by 

2 

"^hi = 3^»(3Ai + A3 + A4 + Ascose). (13) 

The neutral components of the other two doublets ip' and (f" mix into two complex neutral 
states and their masses are given by 

"^n" = Y (^-^3 - A4 - 4A5 cose ± ^J (A3 + A4)2 + 4A2(1 + 2 sin2 e) - 4(A3 + A4)A5 cos 

(14) 

^In the special case where e = 0, the symmetry of the vacuum is enlarged to S3 even if S3 is not a 
subgroup of A4. The reason is that setting e = effectively enlarges the symmetry of the potential to S4 
(once also SU{2) x U(l) gauge invariance is required), which does have S3 as a subgroup. 
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The charged components of </?', do not mix, their masses being 

m';!!^'^ = - Y (3^4 + 3A5 cos e ± v^A5 sin e) . (15) 

4.2 e ^ 0: The Alignment {v, 0, 0) 

In the chosen basis, the vacuum ahgnments (v, 0, 0) preserves the Z2 subgroup of A^. 
As we did with the vacuum ahgnment that conserved the subgroup, in this case it 
is useful to rewrite the scalar potential by performing the following Z2 conserving basis 
transformation 

$1 ^ $1, 

$2 ^ e— /2$2, (16) 

$1 is even under Z2 and behaves hke the standard Higgs doublet, while $2 and $3 are 
odd. For what concerns the neutral states, the 6x6 mass matrix is diagonal in this basis 
and with some degenerated entries: using a notation similar to the 2DHM, we have 

^hi = 2Aiv2 , ml^ = ml = ^(Aa + A4 - K)vl , 

where the last state corresponds to the GB. The charged scalar mass matrix is also diagonal 
with 

"^Ci = = > ^1+ = , (18) 

where the last state corresponds to the GB. The degeneracy in the mass matrices are 
imposed by the residual Z2 symmetry. Contrary to the previous case the neutral scalar 
mass eigenstates are real and not complex. 

4.3 e = 0: The Ahgnment {v, v, 0) 

This vacuum alignment does not preserve any subgroup of A4 and it holds that v — v^/ -\/2- 
Prom the minimum equations we have that 

/^' = -|^^^(4Ai + A3 + A4 + A5). (19) 
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The scalar and pseudoscalar mass eigenvalues are given by 

|^(4Ai + A3 + A4 + A5), 
-X,vl , (20) 
0. 

For the charged sector we have 

m^, = f (A3 - A4 - As) , = (A4 + A5) m^3 = 0. (21) 

For Ag 7^ the alignment {v,v,0) has the correct number of GBs, while for A5 = 
we have an extra massless pesudoscalar. However in both cases, A5 7^ or A5 = 0, the 
conditions m^^ > and m^^ > can not be simultaneously satisfied. This alignment is 
therefore a saddle point of the A4 scalar potential we are studying. 



< = -f (A3 + A4 + A5), 

< = f (A3+A4-3A5), 



4.4 6 = 0: The Alignment (^^1,^^25 ^3) 

This vacuum alignment, as the previous one, does not preserve any subgroup oi A4. A 
part from the condition e = 0, we recall that in this case there is the further constraint 
A3 + A4 + A5 = and A5 may assume both positive and negative values since we have 

reabsorbed in the A5 sign the case e = n. 

Let us define = vf + V2 + = (1 + + r'^)vf with s = V2/V1 and r = v^/vi 
respectively. The mass matrix for the neutral scalar states presents two null eigenvalues- 
as we expected since the condition A3 + A4 + A5 = enlarges the potential symmetry- and 
a massive one 

mhi^ = 2Xivl . (22) 

At the same time the mass matrix for the CP-odd states has one null eigenvalue-the GB 
7r° and two degenerate eigenvalues of mass 

<=< = (A3 + A4)^^ (23) 

Notice that for the special case A5 = we have the constraint A3 = — A4 that implies two 
extra massless pseudoscalars. Finally for the charged scalars we have 

1 

^Ci = = 2^3vl , = (24) 

The total amount of GBs is 5 (7) for the case A5 7^ (A5 = 0), so we have 2 (4) extra 
unwanted GBs: this situation is really problematic. We note that the introduction of 
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terms in the potential that softly break A4 can ameliorate the situation with the Goldstone 
bosons. We will analyse soft A4 breaking terms in more detail in 17 . 



5 CP Non-Conserved Solutions 

In this subsection, we consider vacua that exhibit spontaneous CP violation. This occurs 
if the vev of at least one of the Higgses is inherently complex. A global rotation can always 
absorb one of the three phases of the vevs. 

We note that that the two natural vacua of the previous section {v,v,v) and (f,0,0) 
do not have CP violating analogues as they have only one phase that can be reabsorbed 
to make all vevs real. 

5.1 The Alignment (fie^'^S ^2, 0) 

In this case the third doublet is inert and therefore we are left only with two doublets that 
develop a complex vev and after the redefinition, there is only one phase ui. Taking the 
generic solution (t;ie"^^, ^2, 0) the minimum equations are given by 

vi [coswi[2/i2 + 2\i{vf + + (A3 + \i)v^] + X^v^ cos(e + coi)]] = , 
V2 [(2^2 ^ 2\i{vl + v^) + (A3 + \4)vf + X^vj cos(e + 2uJi)] = , 

(25) 

vi [sinwip/i^ + 2\i{vl + V2) + (A3 + Xijvl] - X^v^sin^e + cui)] = , 
^2^1 sin(e + 2uJi) = . 

The last equation can be solved by e = —2uJi or e = —2uJi + vr. Like in section |4], we can 
absorb the second case by a redefinition of A5. The other three equations reduce to 

vi coswi [2/i2 + 2Xi{vl + w|) + (A3 + Xi)vl + X5VI] = , 

V2 [2/^2 + 2X,{vf + vl) + (A3 + X,)vl + X,vl] = , (26) 

vi sinwi [2/i2 + 2Xi{vf + v^) + (A3 + X^)vl + Agw^] = . 

that are simultaneously solved for f 1 = f 2 = v^j a/2 and 

/i2 = -|^(4Ai + A3 + A4 + A5). (27) 
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The neutral and charged 6x6 mass matrices are quite simple and it is possible having 
analytical expression for the mass eigenvalues. For the neutral sector we have 

= ^^^(-^3 - A4 - As) , ml^ = lvl{4Xi + A3 + A4 + As) , 

<s = 7^^(^3 + A4 - As + 2A5 cos Scui) , ml^ = -X^v^ , (28) 



m 



= |f^^(A3 + A4 - As - 2A5 cos 3a;i) , m^o = , 



and for the charged one we have 



= f(A3 - A4 - As) , = |^(-A4 - As) , ml^^O. (29) 

We see that the mass of the fourth neutral boson selects negative vahics for A5, i.e. the 
second solution e = —2ui + tt. It is interesting to sec that in the CP conserved limits 
cui — )■ (or tt), it is not possible to have both m^^ and m^^ (respectively m^i-) positive, 
but that in the general case, there are points in parameter space where indeed all masses 
are positive. This is in particular clear in the region around cos Soui — 0. 

Finally, as for the CP conserved case, for As = we have two problems: an extra GB 
and we cannot have all positive massive eigenstates. 



5.2 The Alignment (^;ie*'^S ^;2e''^% ^3) 

In this case all the doublets develop a vev i^j 7^ 0, so we may have two physical CP violating 
phases. We have the freedom to take UJ3 — 0. In this case the first derivatives system is 
given by 

vi {cosa;i[2//^ + 2Xi{vf + vl + vj) + (A3 + X4){v^ + vj)] + 

H-Asffg cos(e — coi) + vl cos(e + coi — 2uj2)]} = , 
V2 {cosa;2(2/x2 + 2Xi{vl + + vj) + (A3 + Xi){vl + vl)+ 

+As['i;3Cos(e + a;2) + vl cos{e - 0J2 + 2oji)]} = 0, 
vs {2//2 + 2Xi{vl + vi + vj) + (A3 + X4)ivl + vl)+ 

+A5 [vl cos(e -2uji)+ vl cos(e + 2^2)] } = , (30) 
{sincuip/i^ + 2Xi{vf + vj + f^) + (A3 + X4){vl + vl)] + 

+As[v3 sin(e ~ Ui) — v^ sin(e + cui — 2ci;2)] } = , 
V2 {sina;2(2//2 + 2Xi{vf + vl + v^) + (A3 + X4){vf + vl))+ 

+As[-'i;3 sin(e + 002) + vf sin(e -0J2 + 2uji)]} = , 
V3 [As(— "^1 sin(e — 2uji) + w| sin(e + 2cc;2))] = . 
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The last equation is solved for uj2 = —uji and V2 



Defining = rv and 



vf + V2 + vl = vl, the previous system reduces to the three equations 



v:. 



2(2 + r2) 
2(2 + r2) 

„,2 



(2 + r2 



(4 + 2r2)Ai + (l + r2)(A3 + A4) + 
(4 + 2r2)Ai + (l + r2)(A3 + A4) + 



A. 



COSCJi 

As 
sinwi 



(r^ cos(e — oji) + cos(e + Swi)) 



r sm e 



^1 



sin(e + 3ui)) 



(2 + r2)Ai + A3 + A4 + As cos(e - 2uji 



0. 



(31) 



We can solve the third equation in eq. (31 ) in terms of fi"^ and then the second equation 
in terms of As, giving 



As 



■^[{2 + r^)Xi + A3 + A4 + As cos(e - 2uJi)] , 

(r^ - 1)(A3 + A4) sinwi 
— 1) sin(e — uji) — 2 cos e sin(3a;i) 



(32) 



Then the first equation in eq. ( |3l| ) has two possible solutions, for A4 and e respectively 

(33) 



a) tan e 



A4 — —A3 , 

sin 2uJi + sin Aoji 



cos 2ui — cos Acui 

To test the validity of the solution so far sketched it is necessary to check to be in a 
true minimum of the potential and not to have extra GBs a part from three corresponding 



to the GBs eaten by the gauge bosons. However the relations given in eq. (32 ) and eq. (33 ) 



do not allow to get analytical solutions for the scalar masses in case ii). For this reason 
we will consider only three special limits in this case : r ~ 0, r ~ 1 and r very large. We 
think that these limit situations could be the most interesting ones in the model building 



realizations. Indeed models present in literature 11,12 fall in the third case, r very large. 



5.2.1 Case i) 

In this case the constraints A4 = —A3 puts As to zero and enlarge substantially the sym- 
metries of the potential: we have an accidental 0(3) in the neutral real direction and 
two accidental f/(l)s due to A5 = 0. For this reason the neutral spectrum has 5 massless 
particles, the GB 7r° and 4 other GBs, and only one massive state 



The charged scalars are 



Ci 



m 



C2 



2\iV: 



m. 



C3 



(34) 



(35) 
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The massive states are degenerate as in the CP conserving minima studied in sec. (4.3) 
for the case A5 = 0. 



5.2.2 Case ^^) 

As it is not possible to find analytical solutions, here we will study three special limits of 
case ii. 
• r ~ 



In this case we will neglect terms of order r . From eq. (33) we have that for r ~ 



e ~ -Acui + Ntt , (36) 



thus from eq. ( 32 ) we have 



n"^ = -Xivl - (A3 + A4) 



1 — cos 6uJi 
2 — 4 cos 6ui ' 



1 — 2 cos buJi 

Under these approximations the 6x6 neutral scalar mass matrix gives one null mass 
state, m^o = 0, corresponding to the GB and the following five eigenvalues at leading 
order, given by 

(A3 + A4)/(l - 2cos6wi)t;2 

ml ~ [-2Ai + (4Ai + A3 + A4)(l-cos6u;i)/(l-2cos6a;i)]t;2 (38) 
"^hi ~ - [(A3 + A4) cos Guivl/ (1-2 cos 6ui)] vl , 
ml ~ - [2(A3 + A4) sin^ 3wi/(l - 2 cos Gui)] vl , 

where /[Aj] stays for a linear combination of the adimensional A parameters of the poten- 
tial. The previous neutral spectrum present a lightest state that may be too light to be 
phenomenologically acceptable. Assuming that the A's potential parameters run in the 
'natural' range 0.1 ^ 10 or, somewhat optimistically, 10~^ -f- 10^. For what concerns r 
we are in the limit of ~ 0, so as reference value we may take ~ 10^'^ ^ 10^^. By 
combining these two ranges we find upper bounds 

ml < 200 GeV for A^ ~ 100, ~ 10"^ , 
ml < 25 GeV for A^ ~ 10, ~ 10"^ . (39) 
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Since f[Xi] ~ 100 may be obtained only for very peculiar combinations of the potential 
parameters, the previous estimates indicate that for relative tiny value of r the spectrum 
may present very light neutral states. 

On the contrary, in the charged sector we have the two GBs eaten by the corresponding 
gauge bosons, m^^ = 0, and two complex massive states with masses 

"^Ci ~ "[-^4 + (As + A4Cos6a;i)/(l - 2 cos 6a;i)]f^/2 

(40) 

~ -[2A4 + (A3 + 2A4Cos6a;i)/(l-2cos6a;i)]i;2/2. 

• r ~ 1 

In this limit we may write r ~ 1 + 5 and make an expansion in terms of 6 neglecting 
terms of order 6^. Thus we have 

e ~ 7r/2 - - (JcotSwi + A^TT, (41) 

and then 



;2 



= -(3Xi + Xs + X4)/Svt-6/9{X3 + X4K 
A5 = (5(A3 + A4) cscSwi . (42) 

Under these approximations the 6x6 neutral scalar mass matrix gives the usual null 
mass state, fn'^o, corresponding to the GB and the following five eigenvalues 



<~<~/W0(5>^ 

<~<~-(A3 + A4)/3t;^ (43) 
m^^ ~2(3Ai + A3 + A4)/3t;2 , 

where again /[Aj] stays for a linear combination of the A's potential parameters. A analysis 
similar to the one for the case with r ~ shows that the neutral spectrum may present 
very light states. 

In the charged sector we have the GBs eaten by the gauge bosons and two degenerate 
massive state 

~ ~ -^4/2t;^ . (44) 

• r > 1 

In this case we may perform an expansion in term of 1/r and neglect terms of order 



l/r^. From eq. (33) we have that 

e 2ui + Ntt , (45) 
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and then eq. (32) reduces to 

2 \ 2 

As ~ -(A3 + A4), 

(46) 

Under these approximations we find a massless neutral scalar state, m^o = 0, and the 
other 5 neutral masses are given at leading order by 

- 2\ivl , (47) 

~ ~ (A3 + Xi)vl , 

where once more f[Xi] stays for a linear combination of the A's potential parameters. The 
charged scalar mass matrix is diagonal up to terms of order (9(l/r^) with two massive 
degenerate states 

ml^=ml^=Xsvl/2, (48) 

and the correct number of GBs. 

If we consider now eq. ( 47 ) we see that as for r ~ and r ~ 1 the expressions for m^^ ^ 



say that we may have two very light neutral scalars. Taking as reference values for r the 
range 50 ^ 200 we find 

V^]5GeV(lGeV), (49) 



giving 



ml 2 < 50^ GeV^ for r ~ 50 ^ 



m? 2 < 10^ GeV^ for r ~ 200 , (50) 

where 50(10) GeV may be obtained only for very peculiar combination of the potential 
parameters. In other words we expect that also in the majority of the cases for r in the 
range 50 — 200 we will have mf 2 very light. 

In conclusion, taking into account the SM context and the potential given in eq. ([2]), 



the solution (e*'^^, e"*'^^, r)w^/v2 + with r small, close to 1 or large give rise to very 
light states. Of course this does not mean that these states will be light for any value of 
r but it is a quite strong hint that it is possible that this could be what indeed happens. 
As mentioned before, the addition of soft A4 breaking terms to the potential may help in 
the cases of Goldstone bosons or very light bosons. We will discuss these terms in more 



detail in 17 
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6 Bounds From The Higgs Phenomenology 



In this section we analyse the phenomenology corresponding to the different vacua dis- 
cussed above: unitarity, Z and decays and oblique parameters. In this way we manage 
to constrain the parameter space and, in some cases, to rule out the studied vacuum con- 
figuration. 



6.1 Unitarity 



In this section we account for the tree level unitarity constraints coming from the additional 
scalars present in the theory. We examine the partial wave unitarity for the neutral two- 
particle amplitudes for s ^ M^,M|. We can use the equivalence theorem, so that we 
can compute the amplitudes using only the scalar potential described in eq. Q. In the 
regime of large energies, the only relevant contributions are the quartic couplings in the 



scalar potential 22 25 and then we can write the J = partial wave amplitude oq in 



terms of the tree level amplitude T as 



1 



dcos Q T{s) 



1 



■F[\ 



32n ' ' IGtc 

where F represents a function of the Aj couplings. Using for simplicity the notation 



(51) 



V2 



\ 



(52) 
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we can write the 30 neutral two-particle channels as follows: 



ZaZb h^jA tO 



(53) 



Once written down the full scattering matrix a^, we find a block diagonal structure. The 
first 12 X 12 block concerns the channels 

ZiZi Z2Z2 Z^Z^ III 



Wi Wi , , W3 W3 , 



72 ' ' 72 ' 72 " ' 72 

while the other three 6x6 blocks are related to the channels 



h^Zi , h2Z2 , h^z^ 



wtwb , w^w^ , h'lzb , hb^a , ZaZb , h^h 



b ) 



once we specify the labels (a, h) as (1, 2), (1, 3) and (2, 3). Notice that up this point 
the analysis is completely general and is valid for all the vacua presented. We specify the 
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vacuum configuration, expressing the quartic couplings Aj in terms of tlie masses of tlie 
scalars. Afterwards, putting the constraint that the largest eigenvalues of the scattering 
matrix is in modulus less than 1, we find upper bounds on the scalar masses which we 
use in our numerical analysis. 

6.2 Z And VF^ Decays 

From an experimental point of view gauge bosons decays into scalar particles are detected 
by looking at fermionic channels, such as for example Z — )■ hA — )■ 4/ in the 2HDM, or 
Z decays into partial or total missing energy in a generic new physics scenario. From 
this point of view gauge bosons decays bound the Higgs sector in an extremely model 
dependent way. However since in the SM the Z and the decays into 2 fermions, 4 
fermions or all have been precisely been calculated and measured, we may focus on the 
decays Z, — )■ all. Doing this we overestimate the allowed regions in the parameter 
space, but we have a first and model independent cut arising by the gauge bosons decay. 
Once we will pass to a model dependent analysis the region may only be restricted, not 
enlarged. Furthermore, defining the contribution from new physics as AF, since 

AT%^^ ~ AT%^^ ~ ATf^^ « Tz,w^ , (54) 

we expect the error we commit being quite small. 
From LEP data we have 

^^zw± = r|^± + ATz,w± (55) 



with ATz ~ 0.0023 GeV and ATw± ~ 0.042 GeV |26]. Therefore we may calculate the 
width 

Z — )■ hihj , 



W+ H+hj. 



(56) 



for the different multi Higgs (MH) vacuum configuration studied and select the points that 
satisfy 

r|f^± < Arz,w± . (57) 

Here we have indicated the generic Z — t- hihj referring to our notation introduced in 
section [2] Clearly when CP is conserved the hi have defined CP and only couplings to 
CP odd states are allowed. Of course this is not true for the configuration when CP is 
spontaneously broken. 
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In the vacuum analysis we did we have seen that in few situations we have extra 
massless or very light particles. For those cases the gauge bosons decays put strong 
bounds. For what concerns the Z decays we have 



kz < AFz- 



ky. 1 



mz 
2 \ 3 



<ATz 



I6n Acl 



w 



mi 



mz g 

1677 4cL 

< ATz ^ 

mz g^ 



if both particles hi and hj are massless . 
if hj is masslees and < m^. < m| , 
if hi, hj 7^ and < m^. + m\, < m| . 



(58) 

where g is the SU{2) gauge coupling, cy/ the cosine of the Weinberg angle 6w and the 
parameter kz is given by 



(59) 



with U defined in eq. ([6]). 

Similarly for the decays we have 



k 



w 



m 



2 \ 3 



m 



w 



IGvr 4c 

< ATw ^ if hj is massless and m^. < m 



mw r 



2 

W 



(60) 



where, in analogy to the Z decay, the parameter kw is given by 



ct ttT 



5, 



u, 



T 



{a+3)b^ ia+3)c 



(61) 



with S defined in eq. ([T]). 



6.3 Large Mass Higgs Decay 



Electroweak data analysis considering the data from LEP2 27 and Tevatron 28 put an 



upper bound on the mass of the SM Higgs of 194 GeV at 99% CL 26 . In a MH scenario 



this bound may be roughly translated in the upper bound for the lightest scalar mass, 
m/ij. For large values of the SM Higgs mass, mh > 2mw, the main channel decay is 
h — W^W~ and the upper bound is completely model independent. Let us indicate as 
rvi^(194) the branching ratio of the SM Higgs into two at a mass of 194 GeV. 

In a MH model the lightest Higgs boson couples to the gauge bosons with a coupling 
that is 

ghizz = Pgh¥z^ 

ghiww = P ghww J 



(62) 
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with /3 < 1. In our case for example f3 is given by 

fa (cos Ua U'^i + sin Uj U J+3)i ) , (63) 

with fa = Va/vw and Ua the corresponding CP phase. Taking into account that hi is less 
produced then the SM Higgs and that its T^^{mfiJ is reduced with respect to the SM 
one, 

^wwi^H.) - |/3rr^^(m.J < r^^(194) , (64) 
we can roughly constrain the upper bound for masses rrih^ > 194 GeV. 



6.4 Constraints By Oblique Corrections 

The consistence of a MH model has to be checked also by means of the oblique corrections. 



These corrections can be classified 29 -33 by means of three parameters, namely TSU, 



that maybe written in terms of the physical gauge boson vacuum polarizations as 34 
47r 



T 

S 

u 



9 2 2 
„2 2 

167r^^ 



Azz{ml) - Azm _ ^, ^^^^ _ %:^a;,(o)' 



rrir 



Cw-Sw 



'IQtt 



Awwi^w) - Aww{^) 2 Azz{m%) - AzziO) 



m 



w 



rrir 



- s-'^A'^^iQ) -2swCwA'^z{^) 
(65) 



where s^y, cw are sine and cosine of 6w and e is the electric charge. EW precision mea- 
surements severely constrain the possible values of the three parameters T, S and U . In 
the SM assuming m\ > m| we have 

3 



log 



ml 



m 



2 ' 



SM 



127r 
0. 



log 



mt 



m 



2 ' 



(66) 



For a Higgs boson mass of mh = 117 GeV (and in brackets the difference assuming instead 
mh = 300 GeV), the data allow 26 



gexp ^ 0.10±0.10(-0.08) 
rpexp ^ 0.03±0.11(+0.09) 
Ijexp ^ 0.06±0.10(+0.01) . 



(67) 
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The constraints in eq. (67) must be rescaled not only for the different values of the Higgs 
boson mass but also for a different scalar or fermion field content: for example, if we 
assume to have a MH scenario this gives a contribution to the T-parameter and we 
need 

JlNSS _ rpSM _ rpexp ^gg-J 



A detailed analysis on the TSU in a MH model has been presented in 35 , 36 where all the 
details are carefully explained. However the resulting formulae are valid only for scalar 
masses larger or comparable to rriz- Since this is not the case for a generic MH model and 
particularly for the configurations studied so far, where we have a redundant number of 
massless or extremely light particles, we improved their results, getting full formulae valid 
for any value of the scalar masses (see the appendix A for details). 

7 Results 

We have performed a numerical analysis for all vacuum configurations considered, ne- 
glecting the alignment {v, v, 0) since in this case there are tachyonic states. Our aim was 
to find a region in the parameter space where all the Higgs constraints were satisfied for 
each configuration considered. We have analysed the points generated through subsequent 
constraints, from the weaker one to the stronger according to 

• points Y: true minima -all the squared masses positive- (yellow points in the figures); 

• points B: unitarity bound (blue points); 

• points G: Z and decays (green points); 

• points R: TSU parameters (red points). 

The ratios B/Y, G/B, R/G may be used to indicate which is the stronger constraint 
for each allowed minima. For almost each case we have compared the masses of the two 
lightest neutral states -except for the alignment studied in sec. |5.2.1 where we have only 



one massive neutral state- and the mass of the lightest neutral scalar versus the mass of 
the lightest charged one. Then we have plotted the TS oblique parameters for all the 
green points to check that T is the most constrained one -for this reason we have not 
inserted the plots concerning U. 
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On the contrary for the CP breaking ahgnment (f e*'^^ , f e"*"^^ , rv) we have personahzed 
the plots for reasons that will be clear in the following. 

Notice that in all the following discussion, we refer as mi {1712) to the (next-to-the-) 
lightest neutral state and as rrich^ as the lightest charged mass state. 

7.1 CP Conserved Solutions 
7.1.1 The Alignment (v,v,v) 



In sec. 4.1 we have redefined the initial 3 doublets in term of the surviving symmetry 
representation: 1, 1', 1". One combination corresponds to a singlet doublet, that 
behaves like the SM Higgs: it develops a non-vanishing vev, gives rise to a CP even state 
which we call hi and to the three GBs eaten by the gauge bosons. The others two doublets, 
if' and if", are inert. From these informations we may already figure out what we expect 
by the numerical scan: 

1) when rrihj^ is the smallest mass, hi is the lightest state and corresponds to the SM-like 
Higgs. As a result, the usual SM mass upper bound applies. On the contrary as 
long as we do not consider its coupling with the fermions we do not have a model 
independent lower mass bound. This is due to a combined effect of the CP and Z^ 
symmetries: hi is CP even and singlet under Z^, but couplings like Zhi(f'^, Zhiif'"^, 
W'hiip'^ or W~hiip"^ are forbidden because of Z^ and then gauge boson decays 
cannot constrain the lower mass of hi. 

2) When (f'^ if"^) is the lightest state, we do not have an upper bound on this state 
because the couplings (p'^W~^W~ {ip"^W~^W~) is absent. On the contrary we may 
have a lower bound because couplings like Zip'^^if'"^ and W~ip'^(p"^ are allowed. 

Combining the two situations sketched in points 1) and 2), we expect neither lower nor 
upper bounds for the lightest Higgs mass: according to which of the two cases is most 
favored, we may expect a denser vertical region around mi ~ mz/\/2 when the Z decay 



channel closes according to eq. (58) -case 2) more favored- or a denser vertical line around 
mi ~ 194 GeV, if the large Higgs mass decay constrain applies -case 1) more favored. 
Indeed by looking at fig. [l] we see that we may find R (allowed) points for very tiny mi 
masses and up to ~ 500 GeV when the unitarity bound starts to show its effect. However 
by looking at the crowded points in fig. [T] it seems that case 2) is slightly preferred with 
respect to case 1). Finally for the G points -those that satisfy the minimum, unitarity 
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and decays conditions- we have compared the contributions to the obhque parameters T 
and 5* to see which of the two is more constraining. It turns out to be T, while we have 
not reported U because its behavior is very similar to S. 
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Figure 1: CP conserving alignment {v,v,v): the upper panels show the lightest neutral 
mass mi versus the second lightest neutral mass m2 and the lightest charged one rrichi 
respectively. The gray arc delimits the region below which the Z ( W) decay channel 
opens. On the left plot the arc is only of 45° because m2 > rrii. For points below the arc 
the Z (W ) decay may happens. The points allowed stretch in the region close to the border 



because of the conditions of eq. (58). The dashed vertical lines indicates the approximated 
cuts that occur at mi ~ mz/\/2 and mi ~ 194 GeV according to case 2) and case 1) 
respectively as explained in the text. The down panels show the contributions to T and 
S for the G points. The gray dashed lines indicate the experimental values at 3, 2, la 
level -long, normal, short dashing respectively. The T parameter turns out to be the most 
constraining one. 
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7.1.2 The Alignment (v, 0, 0) 

For what concerns the second natural minimum, the Z2 preserving one, things slightly 



change with respect to the surviving case. By sec. 4.2 we know that as for the case 
we have a SM-like doublet, Z2 even, that develops the vev, gives rise to a CP even neutral 
state, hi, and to the GBs eaten by the gauge bosons. However contrary to the Z^ case, 
in the Z2 minima we have 4 Z2 odd states, 2 CP even labelled /i2 3 and 2 CP odd labelled 



/i4 5. Moreover the 2 CP even (odd) are degenerate. As done in sec. 7.1.1 we may sketch 
what we expect from the numerical analysis: 

1) when hi, the Z2 even SM-like Higgs, is the lightest we expect the SM Higgs upper 
bound but no lower bound because the interactions Zhih^^^ are forbidden by the Z2 
symmetry; 

2) when the two lightest are the Z2 odd degenerate states h2.3 -CP even- or h^^^ -CP 
odd- we expect no upper bound. Moreover since they are degenerate we do not 
expect lower bound too. On the contrary we expect that Z and W decays constrain 
the third lightest neutral Higgs mass and that of the charged ones. 

By looking at fig. [2]we see that indeed we have a large amount of points for which mi = 1712 
for values from up to 700 GeV, thus refiecting case 2). Then the points corresponding 
to case 1) have a sharp cut at mi = 194 GeV, that rejects many blue points, i.e. those 
satisfying the unitarity constrain but not the decays one. We have reported also mi 
versus m^, to check that indeed, when mi — )■ 0, is bounded by m^ as we expected. Our 
intuitions are also confirmed by the plot mi — rrichi- As for the Z3 preserving case the 
most constraining oblique parameter is T. 

7.1.3 The Alignment {vi, V2, v^) with e = 0, A3 + A4 + A5 = 

In this case we do not have any surviving symmetry which forbid some couplings. However 



from sec. 4.3 we know that the conditions e = 0, A3 + A4 + A5 = give rise to two extra 
massless CP even particles. Therefore we expect that 

1) when the lightest massive state is CP odd, then its mass is bounded by the Z decay 



through eq. (58); 



2) when the lightest massive state is CP even, then its mass could reach smaller values 
since the Z decay bound would constrain the combination of its mass with the 
lightest CP odd state mass. 
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Figure 2: CP conserving alignment (f,0,0): the upper panels show mi versus m2 (on the 
left) and third lightest (on the right). For the latter we reported only the R points. The 
central panel shows mi versus m^hi ■ The gray arc delimits the region below which the Z 
(W ) decay channel opens while the second dashed vertical one the SM-Higgs mass upper 
bound at 194 GeV. The first dashed vertical line at mi = mz/V2 is reported to help a 
comparison with the preserving case. On the first two plots the arc is only of 45 degrees 
because m2,3 > mi. The down panels show the contributions to T and S for the G points. 
The T parameter turns out to be the most constraining one. 
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Moreover in both cases we expect the mass of the hghtest charged scalar bounded by W 



decay, according to eq. (60), due to its couphng with W and the massless particles. 



By fig. [3] we see that it seems that case 2) happens very rarely because the cut at 
nil ~ is in evidence. As for the and Z2 preserving minima the T parameter is the 
most constraining one. 
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Figure 3: CP conserving alignment {vi,V2,V3): the upper panels show mi versus m2 and 
fnchi respectively. The dashed lines at mi = mz (vertical) and mchi = i^w (horizontal) 
delimit the region below which the Z and W decay channels open respectively. The allowed 



points concentrate close to the borders according to eqs. 58_ 60_. The down panels show the 
contributions to T and S for the G points. The T parameter turns out to be the most 
constraining one. 



7.2 CP Non-Conserved Solutions 
7.2.1 The Alignment (ve^'^S 0) 



As for the vacuum alignment {vi,V2,v^) commented in sec 7.1.3 the alignment (fe*'^i,f,0) 



does not preserve any subgroup. Moreover since even CP is broken any symmetry 
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cannot help us in sketching the behavior we expect. In general any state, having a CP 
even and a CP odd component, may couple to Z and to another neutral state. However 
we expect limit situations in which for example CP is almost conserved and the 2 lightest 
states have almost the same CP parity. Thus for those cases we do not expect any lower 
bound on m\ and m2. On the contrary the coupling between the W with the lightest 
neutral and the lightest charged scalars does not go to zero when CP is almost restored. 
Then we expect that the quantity + m^^^^ is bounded by the W decay (fig. |4]). For 
what concerns the upper bound on the lighetst neutral mass state we do not expect any 
clear cut because we may not identify a SM-like Higgs. 
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Figure 4: CP no conserving alignment (fe*'^^, f , 0): as in the previous figure the upper 
panels show mi versus m2 and rrichi respectively. In the plot on the right, the effect of the 
W decay constraint on m\ + m^^^ is clear by looking at the B points. The down panels 
show the contributions to T and S for the G points. The T parameter turns out to be the 
most constraining one. 
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7.2.2 The Alignment (ve*'^^, i>e~*'^% rv) case i) 



In sec. 5.2.1 we have seen that the alignment (f e*'^^, fe"*'^^, rf ) with the constrains A5 = 0, 
A4 = —A3, gives rise to 4 extra GBs and only to one neutral state. The simplicity of 
the analytical expressions for the three no vanishing masses ensures that the boundness 
constrain Ai > in addition to A3 > give positive masses. Thus in this case the Y points 
are superfluous. As in the previous cases, we expect the B points to be similar to the Y 
ones, because we choose our parameters centered in 1 in order not to have problems with 
unitarity. In conclusion, for this case only the G and R points are interesting. Moreover 
we expect that the most stringent bound is given by the decay constrains and not by 
TSU: massless particles give a small contribution to the oblique parameters and due to 
the limited number of new particles (2 charged degenerate scalars) TSU should not deviate 
too much by the SM values. Indeed in fig. [5] it is shown that the oblique parameters at 
3 a level do not constrain at all the G points. For this reason we reported only the R 
points in the upper panel of fig. [5j By looking at the plot mi — rrichi in fig. [5] we see that 
with respect to the minima so far analyzed we have much less points and that as expected 
there are cuts in correspondence of mz and mw- 

In conclusion, the solutions for the alignment (fe*^^, f e"*"^^, rf ) with A5 = 0, A4 = —A3 
are not easy to find, but the Higgs phenomenology does not completely rule out this 
vacuum configuration. We could introduce a weight to estimate how much a solution is 
stable or fine-tuned but this goes over the purposes of this work. We expect that this 
situation with 4 extra massless particles could be very problematic when considering the 



model dependent constraints 17 



7.2.3 (t^e^'^S ve-*'^i, ri;) case ii) 



In the analytical discussion done in sec. 5.2.2 we have seen that at least in the special 



limit r ~ (r ~ 1 and r >> 1) we expect the presence of one (two) very light particles. 
From all the numerical scans we performed we found out that solutions for the vacuum 
alignment (f e*"^^ , f e"*"^^ , rf ) with the constraints of case ii) are very difficult to be found. 
Moreover from fig. |6] we see that for any value of r the two lightest states are always very 
light, thus confirming our rough analytical approximations. Indeed both mi and m2 are 
lighter then we expected -especially m2 for r ~ 0- thus indicating that some cancellations 
have to occur to give all the masses greater then 0. This supports the difficulty to find 
solutions, difficulty that cannot to be ascribed to any constrain we imposed, because even 
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Figure 5: CP no conserving alignment {ve'^'^'^, ve~^'^^,rv) case i): the upper panel show mi 
versus rrichi- Only the R points are reported. The down panels show the contributions to 
T and S for the G points. For this specific case the TSU oblique parameter constrain is 
irrelevant compared to the decay one. 
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Figure 6: CP no conserving alignment (f e*"^^, f e"""^^, rf ), case ii): the panels show mi (on 
the left) and m2 (on the right) versus r. The number of points is small, but the interesting 
information is the order of magnitude of the masses. 
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in presence of 4 additional GBs as in sec. 7.2.2 we found out a significant larger number 
of solutions. 

The presence of a single R point in fig. |6] is not statistically relevant, but more in- 
teresting is the order of magnitude of mi^2'- even in case ii) we expect that the align- 
ment (fe*'^^,we~*'^^,rf) may present serious problems once we add model dependent con- 
straints [TtI. 



8 Conclusions 

Flavour models based on non-Abelian discrete symmetries under which the SM scalar 
doublet (and its replicants) transforms non trivially are quite appealing for many reasons. 
First of all there are no new physics scales, since the flavour and the EW symmetries are 
simultaneously broken. Furthermore this kind of models are typically more minimal with 
respect to the ones in which the flavour scale is higher than the EW one: in particular the 
vacuum configuration is simpler and the number of parameters is lower. We then expect an 
high predictive power and clear phenomenological signatures in processes involving both 
fermions and scalars. 

Due to the restricted number of parameters and the abundance of sensitive observables 
in these models, there are many constraints to analyze: the most stringent ones arise by 



FCNC and LFV processes 17 but even Higgs phenomenology put several constraints on 
this class of models. The impact of the symmetry breaking in cosmology has been studied 
in 
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In this paper we focussed on the discrete group, but this analysis can be safely 
generalized for any non-Abelian discrete symmetry. We consider three copies of the SM 
Higgs fields, that transform as a triplet of A4. This setting has already been chosen in 



several papers 10 13 due to the simple vacuum alignment mechanism. 

We have considered all the possible vacuum configurations allowed by the A4 x SM 
scalar potential. These configurations can either conserve or violate CP. For all of them 
we have considered only model independent constraints, related to the Higgs-gauge boson 
Lagrangian, and postponing the model dependent analysis to an accompanying paper [17]. 
The first model independent constraint comes from the partial wave unitarity for the 
neutral two-particle amplitudes, which puts upper bounds on the scalar masses. Then we 
have explained how the light scalar mass region can be constrained considering the gauge 
boson decays. Moreover we have seen how to put an upper bound on the lightest neutral 
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state mass considering the Higgs decay channel h — )■ W^W~ . Finally the most stringent 
bounds arise by the oblique parameters TSU . 

We have shown that the Higgs-gauge boson model independent analysis can be used to 
study the parameter space of the difference vacuum configurations. Among the possible 
solutions which minimize the scalar potential, only one is ruled out due to the presence of 
tachyonic states. Furthermore, some other configurations may be obtained only by tuning 
the potential parameters, giving rise to scalar spectrums characterized by very light or even 
massless particles. Finally, for the remaining ones, we find that they may share common 
features and this increases the difficulty in discriminating among them. Nevertheless, 
the model independent approach restricts in a non trivial way the parameter space. In 
conclusion, we underline that more constraining results can be found considering specific 
realizations which adopt the different vacuum configurations: we present this analysis 
in frf 



Note Added In Proof 



While completing this paper we received ref. 38 , where the scalar potential with three 



copies of the SM Higgs doublet transforming as a triplet of is also studied. We stress the 



differences between this work an ours. Firstly, in 38 , it is assumed that no new CP phases 



appear in the Higgs vevs, while we take this important possibility into account. Secondly, 



ref. 38 discusses three interesting, but rather arbitrary vacua, where our analysis exhausts 
all possible vacua configurations. Lastly, a complete phenomenological study is missing 
in |38 . 
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Appendix A: Analytical Formulae for TSU Parameters 

to our notations 



In this Appendix we provide a sort of translator from the papers 30 , 33 
and furnish the formulae we have used when different from their. 

Reminding their notation we are in the case in which = 3 and nn,nc = so we do 
not have the matrices T and TZ. Then we have 



u 


— > 


S 


ReVfci 




Uki , 


ImVfci 




Uk+3i 









(A.l) 

Moreover they put the GBs as first mass eigenstates while we put them as the last ones 
and contrary to them we use the standard definition for the photon. 
We have rewritten they expression for 



A{I,J,Q)-AiI,J,0) ^ J '^^ 
Q I QFjQ) 

[ Q 



487r2 



for / ^ and/or J ^ , 

logQ for / = J = since A{0, 0, 0) = . 

(A.2) 



For the first row of eq. (A.2) we have used 

dA{I,J,Q) 



A{I,J,Q)^A{I,J,0) + Q- 



dQ 



A{I,J,0) + QdA{I, J) 



(A.3) 



with 



288(/ - J)%^ 



[/3 + 9 J/2 + 6(/ - 3 J) log(/)/2 - 9 J2/ - J3 + 6(3/ - J) J2 log( J)] 



dA(I, J) = < 



1 



2887r2 
1 , 



L487r 



{l + 6Log[I]) 

2vl + l0gW) 



for I,J y^O,I ^ J 

for J = , 
for / = J . 



(A.4) 

The function A{I, J, Q) enters only in the loops in which a gauge boson and a scalar run, 
so we have always J = Q when computing the quantity 

A{I,J,Q)-A{I,J,0) 



Q 



dA{I, J) . 



(A.5) 
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As a result, for this function, it does not make sense considering the case I — J — being 
J — Q = my the gauge boson mass. We found 



1 



dA{I, Q) = < 



8(/ - g)37r2 
~ for / = , 



[Q (-P + 2Q log(7)7 - 2Q \og{Q)I + Q^)] for 7 ^ g, 7 ^ , 







for 7 = g . 



(A.6) 



31 



References 

[1] G. L. Fogli, E. Lisi, A. Marrone, A. Palazzo, and A. M. Rotunno, Phys. Rev. Lett. 
101 (2008) 141801 [arXiv: 0806.2649]. 

[2] T. Schwetz, M. A. Tortola, and J. W. F. Valle, New J. Phys. 10 (2008) 113011 [arXiv: 
0808.2016]. 

[3] M. Maltoni and T. Schwetz, PoS IDM20 08 (2008) 072 [arXiv: 0812.3161]. 

[4] G. L. Fogh, E. Lisi, A. Marrone, A. Palazzo, and A. M. Rotunno, in the proceedings 
of 4th International Workshop on Neutrino Oscillations in Venice: Ten Years after the 
Neutrino Oscillations, Venice, Italy, 15-18 Apr 2008, pp. 21-28 [arXiv: 0809.2936]. 

[5] M. C. Gonzalez-Garcia, M. Maltoni, and J. Salvado, JHEP 04 (2010) 056 [arXiv: 
1001.4524]. 

[6] A. Gando et al, arXiv: 1009.4771. 

[7] P. F. Harrison, D. H. Perkins, and W. G. Scott, Phys. Lett. B530 (2002) 167 [arXiv: 
hep-ph/0202074]. 

[8] Z.-z. Xing, Phys. Lett. B533 (2002) 85-93 [arXiv: hep-ph/0204049]. 

[9] G. AltareUi and F. Feruglio, Rev. Mod. Phys. 82 (2010) 2701-2729 [arXiv: 1002.0211]. 

[10] E. Ma and G. Rajasekaran, Phys. Rev. D64 (2001) 113012 [arXiv: hep-ph/0106291]. 

[11] L. Lavoura and H. Kuhbock, Eur. Phys. J. C55 (2008) 303-308 [arXiv: 0711.0670]. 

[12] S. Morisi and E. Peinado, Phys. Rev. D80 (2009) 113011 [arXiv: 0910.4389]. 

[13] E. Ma, Phys. Rev. D82 (2010) 037301 [arXiv: 1006.3524]. 

[14] E. Ma, Mod. Phys. Lett. A25 (2010) 2215-2221 [arXiv: 0908.3165]. 

[15] M. Hirsch, S. Morisi, E. Peinado, and J. W. F. Valle, (2010) [arXiv: 1007.0871]. 

[16] D. Meloni, S. Morisi, and E. Peinado, (2010) [arXiv: 1011.1371]. 

[17] R. de Adelhart Toorop, F. Bazzocchi, L. Merlo, and A. Paris, arxiv: 1012.2091. 



32 



[18] M. Hamermesh, Group Theory and Its Application to Physical Problems. Addison- 
Wesley Publishing Company, Inc., Reading, Massachusetts (1962). 

[19] J. F. Cornwell, Group Theory in Physics: An Introduction. Academic Press, San 
Diego, CA (1997). 

[20] G. AltareUi and F. Ferugho, Nucl. Phys. B741 (2006) 215-235 [arXiv: hep- 
ph/0512103]. 

[21] J. F. Gunion and H. E. Haber, Phys. Rev. D72 (2005) 095002 [arXiv: hep- 
ph/0506227]. 

[22] M. J. G. Veltman, Acta Phys. Polon. B8 (1977) 475. 

[23] B. W. Lee, C. Quigg, and H. B. Thacker, Phys. Rev. D16 (1977) 1519. 

[24] B. W. Lee, G. Quigg, and H. B. Thacker, Phys. Rev. Lett. 38 (1977) 883-885. 

[25] R. Casalbuoni, D. Dominici, F. Ferugho, and R. Gatto, Nucl. Phys. B299 (1988) 117. 

[26] K. Nakamura et al, J. Phys. G 37 (2010) 075021. 

[27] R. Barate et ai, Phys. Lett. B565 (2003) 61-75 [arXiv: hep-ex/0306033]. 

[28] GDF and D. Collaborations, arXiv: 0903.4001. 

[29] B. Holdom and J. Terning, Phys. Lett. B247 (1990) 88-92. 

[30] M. E. Peskin and T. Takeuchi, Phys. Rev. Lett. 65 (1990) 964-967. 

[31] M. Golden and L. Randall, Nucl. Phys. B361 (1991) 3-23. 

[32] A. Dobado, D. Espriu, and M. J. Herrero, Phys. Lett. B255 (1991) 405-414. 

[33] M. E. Peskin and T. Takeuchi, Phys. Rev. D46 (1992) 381-409. 

[34] 1. Maksymyk, C. P. Burgess, and D. London, Phys. Rev. D50 (1994) 529-535 [arXiv: 
hep-ph/9306267]. 

[35] W. Grimus, L. Lavoura, O. M. Ogreid, and P. Osland, Nucl. Phys. B801 (2008) 81-96 
[arXiv: 0802.4353]. 



33 



[36] W. Grimus, L. Lavoura, O. M. Ogreid, and P. Osland, J. Phys. G35 (2008) 075001 
[arXiv: 0711.4022]. 

[37] M. Y. Khlopov, Cosmoparticle physics, World Scientific (1999). 

[38] A. C. B. Machado, J. C. Montero, and V. Pleitez, (2010) [arXiv: 1011.5855]. 



34 



